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Using the equivalence between the quantum-mechanical probability density in a quantum 
billiard and the Poynting vector in the corresponding microwave system, current distributions 
were studied in a quantum dot like cavity, as well as in a Robnik billiard with A = 0.4, and 
an introduced ferrite cylinder. Spatial auto correlation functions for currents and vorticity 
were studied and compared with predictions from the random-superposition-of-plane-waves 
hypothesis. In addition different types of vortex neighbour spacing distributions were deter- 
mined and compared with theory. 

§1. Introduction 

The large majority of wave functions of chaotic billiards is chaotic, i. e. at any 
point in the system, not too far from the wall, the wave function may be well de- 
scribed by a random-superposition of plane- waves (RSPW)P 

V>(r) = J2 a n^P iknf , (1) 

n 

where modulus k = \k n \ of the incoming wave is fixed, but directions k n /k and am- 
plitudes are considered as random. As an immediate consequence the wave function 
amplitudes are Gaussian distributed, or, equivalently, their squares p = \ip\ 2 are 
Porter-Thomas distributed, 



where A is the billiard area. For the spatial correlation function of the wave function 
amplitudes one obtains a Bessel function, 



C{n,r 2 )= "\)[^;'' =J {kr\ (3) 



(^ *(7T)y>(r 2 )) 

mr)) 

where r = |r*x — t^|. The brackets denote an average over all positions. All these 
features have been demonstrated by McDonald and Kaufman in their influential 
work on stadium wave functionsP^ Very recently Urbina and RichteJ^ succeeded 
in giving a quantum mechanical justification of the approach, generalizing ideas of 
Hortikar, Srednickif^ and Gornyi, Mirlin.* 1 It is impossible to mention all papers 
which have been published hitherto on the subject. The RSPW approach is not 
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restricted to quantum mechanics. This is why experiments using classical waves 
have played an important role, since for a long time they were the only ones with the 
ability to look into the system. Very recently techniques have been developed which 
yield comparable information for electron flow patterns in mesoscopic structures.^ 
The state of the art of the experiments with classical waves up to the year 1999 is 
presented in reference EJ) . 

Most of the experiments with classical waves have been performed in microwave 
resonator^ J^'^ and vibrating solids In one work light propagation through a 
wave guide with distorted cross-section was studied.^ In all cases the predictions 
of the RSPW approach could be verified. It should be noted that in the general case 
there is no one-to-one correspondence to quantum mechanics, thus demonstrating 
the universality of the approach. Therefor similar ideas have been developed in- 
dependently in the context of room acoustics™ Quasi-two-dimensional microwave 
resonators constitute one prominent exception where the equivalence to quantum 
mechanics is complete, including the boundary conditions. This is no longer true 
in three-dimensional resonators. But even here the approach remains valid R One 
only has to superimpose plane electromagnetic waves with the consequence that ex- 
pression (J3J) for the spatial autocorrelation function has to be modified.® Of course 
there are limitations. Since the RSPW approach completely ignores boundary con- 
ditions, deviations are expected and found in regions close to the boundary,^ ^HJ 
or if the wavelengths are not small compared to the system size. 

The approach definitely cannot be applied to wave functions which are scarred 
along periodic orbits or show regular patterns associated with bouncing balls.'™' It 
was shown already by McDonald and Kaufman that for such wave functions the 
wave function amplitudes are not Gaussian distributed. 3 ' If the billiards are open, 
or if time-reversal symmetry is broken, the wave function are complex, and currents 
are present. In microwave experiments wall absorption is another source of currents. 
The quantum-mechanical probability density is given by 



j(r) = ImO*VV0- (4) 

In quasi-two-dimensional electromagnetic cavities there is a one-to-one correspon- 
dence to the Poynting vector making an experimental determination of j(r) feasible 
as wellPS 

The consequences of the RSPW approach for the distribution of currents have 
been studied in particular by Berggren and coworker in a series of papers 
In open systems there are no longer nodal lines but nodal points, or vortices, since for 
the wave function to be zero both real and imaginary part have to be zero at the same 
time. Two-point correlation functions of vortices have been given independently by 
Berry and Dennis^ and by Saichev et al. 2 -' Nearest neighbour distributions of 
vortices have been studied in Ref. I27|). The theoretical predictions have been tested 
experimentally in two microwave experiments including a direct visualization 
of persistent currents well-known from mesoscopic physics. In the present paper a 
number of additional microwave tests of the RSPW hypothesis are presented with 
special emphasis on spatial auto correlation functions of currents and vorticities. To 
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the best of our knowledge such quantities have never been studied before, neither 
theoretically nor experimentally. After a short recapitulation of the experimental 
technique in section [21 analytic expressions for a number of autocorrelation function 
are derived in section |31 Various comparisons between experiment and theory are 
presented in section |3| 



§2. Experiment 



In the description of the experiment we can be short, since all details are available 
from our previous publications for the Robnikp^'^S and the quantum dot billiardP^ 
The quantum dot billiard is an open billiard of rectangular shape with rounded 
corners, and an entrance and an exit wave guide attached at opposite sides (see 
Fig. The second one belongs to the family of the Pascal limacon, or Robnik 
billiards. It can be obtained by a complex mapping of the unit circle by means of 
the function u = z + Xz 1 . For the parameter A = 0.4, used in the experiment, 
the classical phase space is completely chaotic,^ apart from, perhaps, tiny regular 
fractions!*^ A ferrite ring has originally been introduced to break time reversal 
symmetry thus giving rise to persistent currents. Time reversal symmetry is broken 
only for a small frequency region (approximately 3 to 5 GHz) but the absorption 
of the ferrite is present for the whole frequency range used for the data analysis (3 
to 10 GHz). In the present context it is not of relevance whether the origin of the 
currents is break of time-reversal symmetry or absorption. 

The linear dimensions of both cavities is of the order of 25 cm (see the scale in 
Fig- and their height is h=8 mm. For frequencies below 18.75 GHz there is a 
complete equivalence between the quantum mechanical wave function amplitudes ip 
and the electric field E z , where the quantum-mechanical eigenenergy E corresponds 
to the square of the wavenumber k 2 . Measuring the reflection amplitude at one 
antenna, or the transmission amplitude between two antennas, the complete scat- 
tering matrix can be obtained, which for isolated resonances reduces to a billiard 
Breit-Wigner function^ 

IV 



II 



n 



Both eigenfunctions and eigenenergies are slightly modified by the presence of 
the antenna which has been neglected in Eq. (J5J (for an introductory presentation 
see chapter 6 of Ref. EJ). From a transmission measurement thus the wave function 
can be obtained including the sign, whereas a reflection measurement only yields the 
modulus. 

For quasi-two-dimensional systems the Poynting vector S = c/(4nr)E x H reduces 
tcP 

S=^Im[E*(r)VE z (r)], (6) 

illustrating the one-to-one correspondence to the quantum mechanical probability 
density (4) stated above. Typical results for the two billiard systems under study are 
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Fig. 1. Plot of wave function amplitude square l^l 2 , probability density j = Im(tf'*V>f r ), and 
vorticity u = (V a; <Pij)(Vj / !i'i) — (V v ^r)(V x^i) f° r the quantum dot billiard (upper row) at 
v= 8.33 GHz, and the Robnik billiard with a ferrite ring insert (lower row) at v= 7.03 GHz. 
The position of the antenna is indicated by a cross ( x ) . 

In the |^| 2 plot (left column) the intensity is converted into a gray scale, where black corresponds 
to zero intensity. In the vorticity plot (right column) black corresponds to a large positive, and 
white to a large negative sense of rotation. 



shown in Fig. ^ Further examples can be found in our previous publications 
In the right column l^l 2 is plotted in a gray scale for a typical frequency. The middle 
column shows the corresponding flow pattern as obtained from Eq. (jSJ), where the 
arrows reflect the Poynting vectors at the respective points. The right column shows 
a plot of the vorticities, or vortex strengths. The vorticity is, up to the factor 1/2, 
just the curl of the currenlP^ and reduces for two-dimensional systems to 

oj = (V^XV^V/) " (Vy^flXV^i), (7) 

where tpR, are real and imaginary part of the wave function. A plot of the vorticity 
is particularly useful to make the vortex pattern visible as is evident from Fig. ^ 

§3. Theory 

It follows immediately from the RSPW hypothesis, as a consequence of the cen- 
tral limit theorem, that the ip{r) can be treated as Gaussian random variables. They 
obey the well-known property that all higher moments can be expressed in terms of 
the second moment. Thus all distributions of interest can be calculatedJEHEl We 
do not follow this route, however, mainly for pedagogical reasons, but start directly 
from Eq. Q to calculate current and vorticity auto correlation functions. 

Writing k n = fc(cos ip n , sin cp n ), we obtain for the derivatives of the wave function 
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( "' ik^a n sin ip n (8) 

n 



<9y 

Using Eq. (J3J) , it follows for the x component of the current 

jx{r) = k^2<a m (cosip n + cos ip m )e- i(k "- km ^. (9) 

In calculating the autocorrelation functions 

C jx (f{,ri) - (Un)Ur2)) (10) 
we use the assumption that the a n are uncorrelated 

{a* n a m ) = (|a n | 2 ) 5 nm . (11) 

It follows from Eqs. @ and JTUJ) 

C,-.(rl,^) ~ A: 2 |a n | 2 |a m | 2 (cosv9 n + cosp m ) 2 e^-^>(*-*»A 

\n,m I 

' cos 2 ip n e~ i% A le Ck A + (cos p n e~ ik A (cos ¥> n e^\ , (12) 



where r = r\ — r 2 . All averages can be expressed in terms of Bessel functions with 
the result 

CiJ?) = (Mkr){J {kr) - cos 2<pJ 2 (kr)} + 2 cos 2 (p^kr)} 2 ) , (13) 

where we have written Cj x (r) instead of Cj x (f{,f 2 ) to indicate that the autocorrela- 
tion function depends on r = \r\ — r 2 \ exclusively. The normalization Cj x (0) = 1 was 
applied. It only remains to perform the average over ip, the angle between vector r 
and the x axis. The averaging gives 

C jx (r) = [Hkr)f + [J l {kr)f. (14) 

For the autocorrelation function oij y (r) the same expression is obtained. Instead 
of averaging over <p, we may alternatively look for two other quantities, namely the 
autocorrelation functions of j\\(r), and j±(r), the current components parallel and 
perpendicular to r, respectively. Inserting ip = and ip = ^ we obtain from Eq. (fl3|) 



C h {r) = Jo(kr)[J (kr) - J 2 (kr)\ + 2J 1 (Ar) 

C j± (r) = J (kr)[J (kr) + J 2 (kr)] (15) 
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In the same way the spatial autocorrelation function for the vorticity is obtained, 
entering definition (7) with expressions (8) for the derivatives of the wave function. 
The derivation is a step-by-step repetition of the calculation for Cj x (r), and we find 



_ (u{ri)u}(r 2 )) 

= [Jo(At)] 2 - [Jx{kr)f 
§4. Results 



(16) 



It was explained in Ref. I27j) that there is a problem with the experimental 
determination of field and current distributions: the probe antenna moving though 
the billiard gives rise to a leakage current, which is critical close to positions where 
the field amplitudes are large. For the quantum dot billiard in addition there are 
frequencies where the total transmission is zero. In such a case the situation is even 
worse, since now there is exclusively the leakage current from the entrance to the 
probe antenna. All these problematic frequency regions have been omitted from the 
data analysis. 




Fig. 2. Experimental spatial autocorrelation functions of the wave function amplitude for the quan- 
tum dot billiard (left) and the Robnik billiard with ferrite insert (right). The solid lines corre- 
spond to the prediction from the random-superposition of plane- waves approach (see Eq. 0). 

For the quantum dot billiard there is another problem. In the low-frequency 
regime the wave functions are still reminiscent of the rectangle with the checker- 
board patterns typical for such systems. In addition there are frequencies showing 
strong scarring associated with bouncing balls and classical trajectories through the 
system. (The relation between scarring and transport was our original motivation to 
study this systempS) triggered by observations in correspondingly shaped quantum 
dots). This is why for the quantum dot billiard only frequencies above 4.2 GHz were 
considered, where the system behaves chaotically. 

All correlation functions discussed in section |21 depend on the parameter kr ex- 
clusively. Therefore it is possible to superimpose the results from different frequencies 
by an appropriate rescaling to improve statistics. We start with the presentation of 
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Fig. 3. Experimental spatial autocorrelation functions of j x {r), j y (r), j\\{r), and j±(r) (from top 
to bottom) for the quantum dot billiard (left column) and the Robnik billiard with ferrite insert 
(right column). The solid lines correspond to the predictions from the random-superposition of 
plane- waves approach (see Eqs. 1141 to iJTSJ). 



our results for the spatial autocorrelation function of the wave function amplitudes 
(see Fig.|2J). A perfect agreement is found for both systems between the experimental 
results and the prediction from the RSPW hypothesis. This may be considered as a 
check for the validity of the approach in the selected frequency regimes. 

In Fig. |2J the results of the different current autocorrelation functions introduced 
in section |3J are shown for the Robnik billiard with ferrite insert. The corresponding 
figures for the quantum dot billiard have been omitted, since the results for the 
systems are more or less identical. Though there are deviations in detail, the overall 
qualitative agreement is very good. In particular the qualitatively different behaviour 
for the various types of current autocorrelation functions is reproduced correctly. 
For the vorticity autocorrelation function, shown in Fig. 01 the agreement between 
theory and experiment is nearly perfect. One only can speculate why this is the 
case: To determine the current, one needs the product of the wave function with its 
derivative, whereas the vorticity is obtained from the product of two derivatives. As a 
consequence, in the latter case all current offsets slowly varying with the position are 
eliminated. Therefore the vorticity is less sensitive on the mentioned experimental 
imperfections, as it seems. 

In Fig. [51 vortex pair correlation functions for both systems are presented. Pre- 
liminary results for the Robnik billiard have already been shown in Ref. I27|). where 
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Fig. 4. Experimental spatial autocorrelation functions of the vorticity u)(r) for the quantum dot 
billiard (left) and the Robnik billiard with ferrite insert (right). The solid lines correspond to 
the predictions from the random-superposition of plane- waves approach (see Eq. flfif )■ 




< tafl 



Fig. 5. Vortex pair correlation function for the quantum dot billiard (left) and the Robnik billiard 
with ferrite insert (right). The solid line corresponds to the theoretical prediction from the 
random-superposition-of-plane waves approach!^ 



a more detailed discussion of this quantity can be found. The expected oscillatory 
behaviour again is reproduced correctly. The deviations between experiment and 
theory at small distances reflect the experimental resolution. The data were taken 
on a grid of 5 mm side length with the consequence that vortices with a distance 
below 10 mm can no longer be separated reliably. 

We now extend this discussion to the nearest neighbour spacing distribution be- 
tween vortices. This quantity was introduced by Saichev et alP^ and was studied by 
the authors in a number of papers. There are different types of spacing distributions 

denoted by P ++ (r),P H (r),P_ + (r),P (r) where the pair of indices denotes the 

sense of rotation of the vortices considered. P++(r) and P (r), as well as P_| (r) 

and P (-(r) should be identical, of course. Only for small system sizes (where the 

RSPW approach fails anyway) there may be deviations due to the presence of the 
boundaryPS Fig. shows our results, for P ++ (r) and P-\ (r). Within the lim- 
its of statistical errors there was no difference to the corresponding distributions 
of P (r) and P \-(r), respectively. The theoretical curve is the result of direct 
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Fig. 6. Nearest neighbour distance distribution of vortices of opposite (P^ (r)), and same sense 

of rotation (P++(r)), respectively, for the quantum dot billiard (left), and the Robnik billiard 
with ferrite insert (right). The solid lines corresponds to the numerical evaluation of the nodes 
in the Berry function. 

numerical evaluation of the nodes calculated via RSPW approach (setPMSB). The 

deviations between experiment and theory for (r) are comparable to that found 

in the papers of the Berggren group P^J^ and reflect at least partly the limitations 
of the Poisson approximation. Another cause, in particular for the small distances, 
is the limited experimental resolution discussed above. For P ++ (r), on the other 
hand, the agreement between experiment and theory is good, again in accordance 
with Refs. I23 |) . I27|) . In addition, the experiment resolution has only a small effect in 
this case, since small distances do not contribute anyway to P++(r) significantly. 

§5. Conclusion 

In this paper a number of consequences of the RSPW hypothesis have been 
presented not been studied hitherto. A qualitatively good agreement between the 
experiment and the theoretical prediction was found for different types of spatial 
current correlation functions. The remaining discrepancies are probably due to ex- 
perimental imperfections caused by leakage currents into the probe antennas. For 
the vorticity autocorrelation function a perfect agreement between experiment and 
theory is found. In addition vortex pair correlation functions, and vortex nearest 
neighbour distance distributions were studied. Again the agreement between exper- 
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iment and theory are good, apart from discrepancies at small distances. 
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